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ABSTRACT 

The two components of radio emission, above and below 86 GHz respectively, from the Galactic center 
source— Sgr A* can be naturally explained by the hybrid of thermal and nonthermal electrons in hot 
accretion flows (e.g., radiatively inefhcient accretion flow, RIAF, e.g.. Yuan et al. 2003). We further apply 
this model to a sample of nearby low-luminosity active galactic nuclei (LLAGNs), which are also believed 
to be powered by the RIAF. We selected the LLAGNs with only compact radio cores according to the 
high-resolution radio observations, and the sources observed with jets or jet-like features are excluded. 
We find that the radio emission of LLAGNs is severely underpredicted by pure RIAF model, and can 
be naturally explained by the RIAF model with a hybrid electron population consisting of both thermal 
and nonthermal particles. Our model can roughly reproduce the observed anti-correlation between the 
mass-corrected radio loudness and Eddington ratio for the LLAGNs in our sample. We further model 
the spectral energy distributions of each source in our sample, and find that roughly all sources can be 
well fitted if a small fraction of the steady state electron energy is ejected into the nonthermal electrons. 
The size of radio emission region of our model is around several thousand gravitational radii, which is 
also roughly consistent with the recent high-resolution VLBI observations for some nearby LLAGNs. 

Subject headings: accretion, accretion disks-black hole physics-galaxies: active-radiation mechanisms: 
nonthermal 



1. INTRODUCTION 

Active galactic nuclei (AGNs) produce enormous lumi- 
nosities in extremely compact volumes. It was found some 
local galaxies show similarities to those of bright AGNs, 
but with relatively weaker broad-li ne emissi on and mul- 
tiwavelength radi ative c ores (e.g.. Ho et al.l [19 97. 2001; 
iFalcke fc MaTkoa [2000t iNaear et all l2002t iMaozi i2007i) . 
These low-luminosity active galactic nuclei (LLAGNs) are 
very common in the local universe, which may be the 
scale d-down luminosity version of bright AGNs (see iHol 
120081 for a recent review). Actually, the AGN may only 
spend a small fraction of its lifetime in the highly lumi- 
nous, QSO-like phase, and much more time is spent in 
a w eakly accreting state (LLA GNs or normal galaxies, 
e.g.. lHopkins &: Hernquis^l2006( ). Both bright AGNs and 
LLAGNs are believed to be powered by the matter falling 
onto a super massive black hole (BH). The observational 
properties of LLAGNs are quite different from those of 
bright QSOs, which are thought to be driven by differ- 
ent accretion modes. The optical/UV bumps observed in 
QSOs can be naturally interpreted as multi-temperature 
blackbody emission from a cold, opticall y thick, and ge- 
omet rically thin standard disk (SSD, Sha kura fc SunvaevI 
Il973() . However, most of LLAGNs lack the evident op- 
tical/UV bump as that of bright AGNs. The possible 
physical reason is that the standard thin accretion disk 
is absent in these LLAGNs. A hot, optically thin, geo- 
metrically thick radiatively inefficient accretion flow model 
has been developed in the last several decades (RIAF, an 
'updated' version of the original advection dominated ac- 
cretion flow model, ADAF; e.g., Ichimaru 1977; Narayan 
&; Yi 1994, 1995; Abramowicz et al. 1995; see Kato et al. 
2008; Narayan & McClintock 2008; and Yuan & Narayan 



2013 for reviews). The RIAF model can successfully ex- 
plain most observational features of nearby LLAGNs (e.g., 
Quataert et al. 1999; Cao & Rawlings 2004; Nemmen et 
al. 2006; Wu et al. 2007; Wu & Gu 2008; Gu & Cao 2009; 
Yuan et al. 2009a; Yuan et al. 2009b; Xu & Cao 2009; Ho 
2009; Cao 2010; Yu et al. 2011; Nemmen et al. 2012; see 
Yuan 2007 and Ho 2008 for recent reviews). 

The AGNs are traditionally divided into radio loud 
(RL) and radio quiet (RQ), and the origin for this di- 
chotomy is still a matter of debate. The distinction be- 
tween the RL and RQ objects is normally based on a 
radio-loudness parameter, which was defined as the ra- 
tio of the monochromatic fiux density at 5 GHz and op- 
tical B band at 4400 A (i?o = ghz/^b)- The loud- 
ness Ro = 10 is usually taken as the division between 
RQ and RL AGNs (p ar ticularly in QSO studies, e.g., 
iKellermann et al.l 11994) . iTerashima k. WilsonI ()2003[ ) in- 
troduced a new definition of the radio loudness parame- 
ter by comparing the 5 GHz radio luminosity to the 2-10 
keV luminosity, _Rx — -^5GHz/-^^2-iokeV: and proposed that 
log i?x — —4.5 would be the barrier separating RL and RQ 
AGNs, which roughly corresponding to Ro = 10. The use 
of X-ray luminosity with respect the optical one should 
largely avoid extinction problems which normally occur in 
the optical band. However, most of the traditionally RQ 
LLAGNs become RL according to the criterion Ro = 10 
or logi?x = —4.5 if using their n uclear eniission at radio 
and o pt ical/X-ray waveband fe.g.. lHoll2002tlPanessa et all 
|2006() . iPanessa et al.l (|2006() redefined the boundary for 
the RL and RQ LLAGNs as logi?x ^ —2.8 based on a 
sample of low-luminosity Seyferts and Fanaroff- Riley type 
I radio galaxies (FR Is). In the last decade, there is great 
progress in estimating the BH mass in both normal galax- 
ies and AGNs (e.g., iGebhardt etall l2000t iKaspi et al.l 
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|2005() . The measurements of BH mass help us to further 
understand the issue of radio loudness. |Hg ()2002[ ) found 
that the loudness para meter increases w ith decreasing of 
the Eddington ratio. Si kora et al.l (|2007t ) investigated the 
radio loudness of a total 199 sources which include broad- 
line radio galaxies (BLRGs), RL QSOs, Seyferts, low- 
ionization nuclear emission- line region galaxies (LINERs) , 
PR Is, and Palomar-Green QSOs (PGQs), and found that 
there are two distinct, approximately parallel tracks on 
radio-loudness— Eddington-ratio plane. They further pro- 
posed a quantitative definition for the radio-loudness pa- 
rameter that is dependent on the Eddington ratio (see 
their equation 5), where the RL sources include BLRGs, 
RL QSOs, PR Is, while RQ sources comprise Seyferts, LIN- 
ERs and PGQs (hereafter we use this criteria to distin- 
guish the RL/RQ). 

The radio emission of RL AGNs is believed to originate 
in the relativistic jets, where many large-scale radio jets 
have been observed directly. Compared with RL AGNs, 
the origin of nuclear radio emission from RQ objects is 
still not well understood, since that RQ AGNs are much 
fainter radio sources, and most of them only show a com- 
pact radio core. The radio variability of RQ QSOs con- 
fined the radio emission to ^0.1 pc fe.g.. iBarvainis et al.l 
l2005f ). and the high-resolution VLBI observations show 
that the region of the radio emiss ion in RQ LLAGNs 



should be 5, 10~ _ to < IQ-^ pc (e.g..lHo k Ulvestad 2001, 



iNagar et ai.il2002t iGiroletti fc Panessall2009D . Several pos 
sibilities have been proposed to explain the radio emission 
in these RQ AGNs. The first explanation is that RQ AGNs 
are scaled-down version o f the RL AGNs where exist a 
weak , small-scale jet (e.g.. iMiller et al.lll993l: iPalcke et all 
Il996f ). A small fraction of RQ objects show core-jet or 
linear structures at pc-scale or su b-pc sc ale, which do 
support the weak-jet scenario (e.g iHo fc U lvestad 2001; 
lUlvestad et"alll2005t ITeipski et al.l 120061) . The accretion- 
jet model is frequently used to model the multiwave- 
length spectral energy distri bution (SEP) of AGNs (par- 
ticularly in LLAGNs, e.g ., lYuan et al.l l2009al: lYu et al.l 
l2011l : [Nemmen et al.ll2012i i. However, this scenario cannot 
ans wer the bimodality of RL/RQ distribution for AGNs 
(e.g. iKellermann et al.|[T99llSikora et al.ll2007^ . The sec- 
ond possibility is that the radio emission of the RQ AGNs 
come from the hot optically thin plasma in the disk winds, 
where the plasma is completely ionized and has a density 
high enough for bremsstrahlung emission that make a sig- 
nifica nt contribution to the obse rved compact radio cor e 
(e.g., iBlundell fc Kuncic 2006; S teenbrugge et al.l 1201 iD . 
In this model the mass-loss rate in the winds are signifi- 
cant and the observed luminosities of radio emission from 
quasars imply that they should accrete at super-Eddington 
rates. However, most of RQ AGNs are accreting at sub- 
Eddington rates, where the radio core emission should not 
dominantly come from the disk winds. The third possi- 
bility is that the radio emission also com e from the ac- 
cretio n flows (e.g., RIAP or disk corona). iLaor fc Behaij 
(|2008[ ) found that the tight correlation between the ra- 
dio and X-ray luminosities in RQ QSOs is similar to that 
of active stars (I/r ~ 10~^Lx)- They, therefore, pro- 
posed that both the radio and X-ray emission from the 
nuclei of RQ AGNs may d ominated by the m agnetic re- 
connection heated corona. iWu fc Cad (l2005h found that 



the radio emission of most nearby LLAGNs is higher than 
that predicted by radiation from the pure thermal elec- 
trons in the RIAF, and their radio emission should have 
other origin (e.g., nonthermal electrons in RIAP or jet). 
The radio spectrum of Sgr A*, a supermassive BH in our 
galaxy, consists of two components, which domin a te be- 
low and above 86 GHz, respectively. lYuan et al] (|2003() 
found that the component above 86 GHz can be well ex- 
plained by the thermal electrons from the RIAP, while 
low-frequency radio spectrum can be explained if there 
exist a s mall fraction of n onthermal electrons in the RIAP 
(see also lOzel et"aIll2000D . 

It is natural that both the thermal and nonthermal elec- 
trons exist in the hot plasma (RIAP or disk corona) , since 
that the turbulence, magnetic reconnection, and weak 
shocks can accelerate electrons and generate a nonther- 
mal tail at high energies in the distribution function of 
thermal electrons. The Sgr A* provide an excellent ob- 
servational evidence for the existence of hybrid of thermal 
and nonthermal electrons in accretion flows (e.g., RIAP), 
since that no evident radio jet was observed up to now 
and, therefore, its radio emission may be dominated by 
other mechanisms. It is interesting to note that there is 
also independent observatio nal evidence for nonther mal 
electrons in accretion flows. iMcConnell et al.| (j2000t ) re- 
ported a high-energy tail in the hard state of the X-ray 
binary, Cygnus X-1, extending from 50 keV to ~5 MeV. 
The data at MeV energies, collected with the COMPTEL 
instrument of the Compton Gamma-Ray Observatory, can 
be explained by a power-law distri bution of nontherma l 
electrons in the RIAP/corona (e.g.. iRomero et all 120101 ). 
In this work, we try to explore whether the nonthermal 
elctrons in RIAP can explain the radio emission of nearby 
LLAGNs, which always have only compact radio cores. 

2. SAMPLE 

To examine whether the radio emission of nearby 
LLAGNs with only compact cores can be explained by 
the nonthermal electrons in the RIAP, we search the liter- 
atures for high-resolution radio and X-ray data to ensure 
that their emission being from the nuclei of sources. In 
the radio band, we only selected the sources that have 
been observed by high- resolution radio telescopes (e.g., 
VLA, MERLIN, VLB A and VLBI) , which are mainly cho- 
sen from IHo fc UlvestadI (|200l[ ). fpilho et al.l ()2006[ ) and 
iNagar et al.l ( 2005 ). where all these works tried to give 
a complete radio imaging survey of all nearby LLAGN s 
given in Palomar spectroscopic survey (jHo et al.l |1997( ). 
Por purpose of our study, we exclude the sources that 
observed with radio jet or even linear radio structures 
where the radio emission may mainly originate in the 
jets/outflows. The 5 GHz nuclear radio luminosities are 
shown in Table (1), where two sources observed at other 
waveband (NGC 1097 at 8.4 GHz and NGC 4736 at 15 
GHz) w as converted to 5 GHz by assuming F^, cx as 
found in lHo et al.l ()200lD for LLAGNs. In the X-ray wave- 
band, we only selected the sources have compact X-ray 
cores, which are observed by high-resolution telescopes of 
Chandra and/or XMM — Newton. Several sources with 
the high-resolution optical da ta from HST are als o in- 
cluded in building t he SEDs ([E racleous et 311120101 and 
references therein). iMaoz et al.i (,2005.) gave lower limits 
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Table 1 
The sample of LLAGNs 



Name 




log Mbh 


log isGHz 


_Lel6scopG 


logL2-10kcV 


Telescope 


^^out 


V 


jrv.cis. 




Mpc 


Mr-, 


crg/s 




crg/s 












(2) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


noi 




14.5 


8.10 


35.98 


VLA 


40.63 


Q 


1.1 X 10""^ 


0.03% 


12 3 


ly vjr V_.' ^ 1 o 1 


13.0 


8.07 


37.06 


VLBA 


39.28 


Q 




10.0% 


4 4 5 


IN V-T V_y O ±4l I 


40.9 


8.60 


37.78 


VLBA 


41.88 


Q 


5.3x 10"*^ 


5% 


4 4 6 




19.7 


7.79 


37.16 


VLBA 


41.41 


Q 


6.0x 10""^ 


3.0% 


4 4 6 




23.4 


8.07 


37.06 


VLBA 


40.74 


Q 


1.1 X 10""^ 


3.0% 


4 4 g 




20.6 


7.44 


36.31 


VLBA 


41.70 


X 


2.1x 10~^ 


4% 


4 4 7 




24.9 


8.42 


36.59 


MERLIN 


39.89 


Q 


4 2x 10"* 


0.3% 


8,8,5 




17.0 


7.85 


36.77 


VLBA 


40.44 


Q 


1.4x10"'^ 


3.0% 


4,4,9 




18.9 


7.55 


35.69 


VLA 


39.27 


X 


5.5x 10"'^ 


5% 


10 10 7 


NGC 3998 


21.6 


8.72 


38.34 


VLBA 


41.79 


X 


2.3x10-3 


3.0% 


4,4,11 


NGC 4138 


17.0 


7.51 


36.13 


VLA 


41.48 


X 


1.5x10-2 


0.1% 


10,10,7 


NGC 4143 


17.0 


8.25 


37.15 


VLBA 


40.04 


C 


5.0x10-"* 


3.0% 


4,4,6 


NGC 4168 


16.8 


7.98 


36.63 


VLBA 


39.32 


X 


3.8x10-* 


4.0% 


4,4,12 


NGC 4203 


9.7 


7.77 


36.70 


VLBA 


39.70 


C 


6.0x10-* 


5.0% 


4,4,13 


NGC 4477 


16.8 


8.00 


35.46 


VLA 


39.60 


X 


4.95x10-* 


0.04% 


10,10,7 


NGC 4565 


9.7 


7.46 


36.26 


VLBA 


39.56 


C 


7.0x10-* 


7.0% 


4,4,6 


NGC 4594 


9.2 


8.46 


37.79 


VLA 


39.95 


C 


3.2x10-* 


15.0% 


14,15,5 


NGC 4639 


16.8 


6.68 


35.57 


VLA 


40.18 


C 


4.5x10-3 


4.0% 


10,10,13 


NGC 4698 


16.8 


7.41 


35.64 


VLA 


38.70 


C 


3.5x10-* 


3.0% 


10,10,5 


NGC 4736 


4.3 


6.98 


35.51 


VLA 


38.77 


C 


5.7x10-* 


3.0% 


1,4,16 



Note. — Column 1: galaxy name; Column 2: galaxy distance; Column 3: BH mass; Column 4: monochromatic power at 5 GHz; Column 5: 
tlie telescope for radio data; Column 6: 2-10 keV X-ray luminosity; Column 7: the telescope for X-ray data, where 'C denote Chandra and 
'X' denote XMM — Newton; Column 8 and 9: model parameters resulting from the SED fits; Column 10: references for the distance, radio 
luminosity and X-ray luminosity respe ctive ly. ^ ^ ^ 

Refere nces: (1) 'Eracleous ct al. ('201^); (2) 'Thean ct al.' (200U); (3) 'Ncmmcn et al.' [20M); (4) 'Nagar ct al.' ('2005'); (51 'Gonzalez-Martin eT^ 
ll2009l'):C6')ribra shima & Wilson (2003); (7) Caopi ct al. (2006); (8) Filho ct al. (2006); (9) Satyapal ct al. (2005); (10) Ho & Ulvcstad (2 00lf)7 

11) | Ptak eTal.1 \2004 ): (12) iPanessa et al.i \200& ): (131 , Ho et al.. (i2001j); (14) ifla U99a); (15) iHummel et al.i ((1984 1: (16) lEracleous et al.l 

200f r 



for the intrinsic AGN optical emission for a sample of 
LINERs, which were constrained from their optical vari- 
abilities. We also include the lower limits of the intrinsic 
AGN optical emission for three so urces (NGC 3998, NGC 
4594 and NGC 4736, [200?). We estimate the BH 

mass from the velocity dispersion cr* of the host bulges for 
L LAGNs in our sa mple using the A/bh — cr* relation given 
bv lTremaine et all ([2002). Most of velocity dispersions are 
selected from the Hypcrlcda database^ with one exception, 
NGC 1097, from .Eracleous et al, (,2010.1 . We find that the 
velocity dispersion of most LLAGNs in our s ample is con- 
sistent with that reported in lHo et"aLl (|2009( ). Our sample 
include 20 LLAGNs (see Table 1). 

3. MODEL 

In this work, w e consider the RIAF model proposed by 
lYuan et al] (|2003D , which can be con sidered as an update d 
version of the orig inal ADAF (e.g., iNaravan fc Yil 119951 ). 
where both outflows and the possible existence of nonther- 
mal electrons are considered. Here, we briefly summarize 
the model as follows. 

The more accurate global structure and dynamics of 
the accretion flow is calculated numerically to obtain 
the ion and electron temperature, density at each ra- 
diu s. In p articular, we employ the approach suggested 
by iManmot o (2000) for calculating the structure of the 
RIAF in general relativistic frame, which allows us to 
calculate the structure of the accretion flow surround- 



ing either a spinning or a nonspinning BH. In this 
work, we calculate the global structure of the accre- 
tion flows surrounding massive Schwarzschild black holes. 
In stead of using a constant accretion rate, we assume 
that the accretion rate is a function of radius, e.g.. 



M = Mont(R/RontY fe.E.. Blandford & Beeelmanl 


1999 


Stone et al.' 


19991 llgumenshchev & Abramowicz 


1999 


Stone & Prinalc 20011: Li & Cao 2009: Yuan & Bull20100 



where i?out is the outer radius of the RIAF and Mout is 
the accretion rate at i?out- The global structure of the 
RIAF can be calculated with proper outer bou ndaries, if 
the p arameters m, a, /?, and S are specified (see lManmotol 
[2000. for more details). All radiation processes (syn- 
chrotron, bremsstrahlung and Compton scattering) are in- 
cluded self-consistently in the calculations of the RIAF 
structure. The parameter m = M/MEdd is dimension- 
less accretion rate, where the Eddington accretion rate is 
defined as MEdd = 1-4 x IO^^Mbh/M© gs"^ For the vis- 
cosity parameter a and magnetic parameter /3 = Pg/Ptot 
[defined as the ratio of the gas to the total pressure (sum 
of gas and magnetic pressure)], we adopt typical values of 
a = 0.3 and (3 — 0.9, which are widely used in RIAF mod- 
els. There is obviously a degeneracy between rriout and 
s when the accretion rate at the innermost region of the 
RIAF is concerned. We adopt s = 0.3 in this work, which 
is well constrained from t he observation of our Galactic 
center black hole, Sgr A* (jYuan et al.ll2003l but see also 



2 http://leda.univ-lyonl.fr 
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lYuan et al] (|2012l) for slightly higher value ofs'~0.4-0.5 
from simulations). We fix the outer boundary at 5000i?g 
{Rg = GMbh/c^ is gravitational radius), which is roughly 
consistent with the prediction of the disk-corona evapo- 
ration model for the LLAGNs with very lo w Eddington 
ratios (e.g., 10~^ < L2-10 kcV /^Edd < lO""^. iLiu fc TaamI 
l2009f ). The most poorly constrained parameter is 6, which 
describes the fraction of the turbulent dis sipation that di- 
rectly heats the electrons in the flow. .Sharma et al.l (|2007t ) 
found that the parameter i5 may be in range of^ 0.01 — 0.3 
based on the simulations, depending on the model details. 
In this work, we adopt (5 = 0.1 because we find that the 
X-ray spectra of most LLAGNs in our sample can be well 
fitted with t his fixed S v a lue. 

Following lYuan et al.l (|2003D . we assume that the in- 
jected energy in the nonthermal electrons is equal to a 
fraction rj of the energy in thermal electrons, where we 
take rj to be independent of the radius. The thermal elec- 
trons are assumed to have relativistic Maxwell-Boltzmann 
distribution, while the nonthermal electrons are assumed 
to be the power-law tail of the thermal electrons, which 
can be described by the parameter p [n{'j) cx 7~^, where 
7 is the electron Lorentz factor]. The number density of 
the nonthermal electrons can be obtained if r] is given. 
lOzel et al] ()2000D proposed that the parameters, p and ry, 
for describing the nonthermal electrons are degenerate in 
the radio waveband in our model. Therefore, we adopt a 
fiducial value oi p = 3.0, and allow the parameter 77 to be 
free. 

After determining the distribution of thermal and non- 
thermal electrons, we calculate their radiation. The syn- 
chrotron emission from bo th thermal and no nthermal elec- 
trons are calculated (e.g.. lYuan et al.l [20031 ) . The Comp- 
tonization of seed photons from both thermal and nonther- 
mal electrons are co nsidered, which was calcula ted by the 
method proposed by ICoppi fc BlandfordI (|1990( ) . In spec- 
tral calculations, the gravitational redshift effect is con- 
sidered, while the relativistic optics near the black hole is 
neglected. More detail s about spectrum calculation can 
be found in I Yuan et al.l (|2003„. and references therein) . 

In summary, our model have only two free parameters, 
rhout and rj in fitting the SED of LLAGNs. 

4. RESULTS 

We present the typical spectrum of the RIAF with a 
population of hybrid electrons in Figure [T] (top panel), 
where the short-dashed, long-dashed and solid lines repre- 
sent the emission from the thermal electrons, nonthermal 
electrons and sum of them respectively, where typical pa- 
rameters Mbh = 1O^M0, TOout — 10^^, and ?] — 1% are 
adopted for given outer radius i?out = 5000i?g. In this 
model, the X-ray emission mainly come from the Comp- 
tonization of seed photons that produced by both ther- 
mal electrons and nonthermal electrons. We find that 
the low-frequency radio emission is dominated by the self- 
absorbed synchrotron emission from the nonthermal elec- 
trons, which is normally 2-3 order of magnitude higher 
than that from thermal electrons in RIAF. Our results 
are similar for other parameters if the typical parameter 

J I ^ several perc ent is adopted as that used in Sgr A* 
Yuan et al.l[2003L[2006h . The high -energy X-ray emission 
dominantly come from the inner region of the accretion 



flows (e.g., within several tens gravitational radii). How- 
ever, we find the region of the radio emission is much 
larger. The dotted lines in Figure [1] (top panel) show 
the spectrum from the nonthermal electrons with outer 
boundary i?out = 3000i?g, 1000i?g and 500i?g respectively 
for above given other parameters. We find that the radia- 
tion region of the low-frequency radio emission at several 
GHz can up to several thousands gravitational radii, where 
the radio emission at > 3000i?g is negligible. The bottom 
panel of Figure [T] show the spectrum of nonthermal elec- 
trons with p = 2.5, r] = 0.5% (sohd line); p = 3.0, 7] = 1% 
(dashed line) and p = 3.5, rj = 5% (dotted line), where 
the other parameters keep the same as above. It can be 
found that p and 77 are degenerate for radio emission (e.g., 
< 10 GHz, see also lOzel et aI1l2000f) . Therefore, it should 
be impossible to determine p and 77 from the observational 
low-frequency radio emission alone. 
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log(v/Hz) 



Fig. 1. — The top panel show the typical spectrum of the RIAF 
with both thermal and nonthermal electrons for Mbh = 1O*M0, 
™.out = 10""^, r] = 1% and p = 3.0 at given iJout = 5000-Rg, where 
the short-dashed, long-dashed and solid lines are the radiation from 
the thermal electrons, nonthermal electrons and sum of them re- 
spectively. The dotted lines represent the radiation of the nonther- 
mal electrons in the RIAF with different outer radius at 3000-Rg, 
1000-Rg, and 500-Rg (from top to bottom) respectively. The bot- 
tom panel show the spectrum from the nonthermal electrons with 
p = 2.5, ri = 0.5% (solid line); p = 3.0, r/ = 1% (dashed line), 
and p = 3.5, t] = 5% (dotted line) respectively, where the other 
parameters are the same as above. 

Our calculations show that the X-ray emission from the 
RIAF is almost proportional to the BH mass for a given rn, 
and it scales as m'^ with q :^ 2 for a giv en BH mass (see also 
iWu k Gag .200fe .Merloni et al.ll2003[ ). In this work, we 
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further explore the possible relation between the radio lu- 
minosity and BH mass/accretion rate for RIAF model with 
both thermal and nonthermal electrons, where the radio 
emission dominantly originate in the nonthermal electrons. 
Figure [5] (top panel) shows the relation between radio lu- 
minosity and BH mass for rhout = 10~^ and rhout = 10"'* 
with 7] = 1% respectively, where we find that the radio 
emission is roughly proportional to Afg^ (solid lines) . The 
bottom panel of Figure [5] shows the relation between the 
radio luminosity and dimensionless accretion rate for given 
BH masses Mbh = 10^ and IO^Mq respectively, where we 
find that the radio emission is roughly proportional to 
with ^ ~ 0.8 for both cases. To explore the radio emis- 
sion for a sample of LLAGNs with different BH masses 
as a whole, we define a mass-corrected radio-loudness as 
-Rm = Rx/M^^, according to the dependence of the radio 
and X-ray emission on the BH mass in our model, where 

= i5GHz/-^2-10kcV- 




6 7 8 9 10 

Iog(MBH/A^o) 




log(m out) 

Fig. 2. — Top panel show the model prediction on the relation be- 
tween -LsGHz and A/bh for two given accretion rates rhout = 10~^ 
and rhout = 10 ~* with r] = 1%, where the solid points represent 
Mbh = 10*^, 10'^, 10*, 10^, and 10^° Mq respectively, and the solid 
lines are their best fittings. The bottom panel present the relation 
between Lsghz and rhout for two given BH masses A/bh = lO'^ 
and A/bh = W^Mq, where the solid points represent rhout=10~^, 
10~^, lO"*, 10"'', and 10~^ respectively, and solid lines are their 
best fittings. 



We present the relation between the mass-corrected 
loudness, i?M, and the Eddington ratio, L2-iokcv/-^Edd! 
for all LLAGNs in Figure [31 where both quantities roughly 
do not depend on the BH mass. It can be found that there 



is an anti-correlation between Rm and -/j2-iokcv/-^Edd for 
LLAGNs. We further demonstrate the theoretical relation 
between Rm and L2-iokcv/-ZjEdd in our model for different 
values of t] from to 10% at given Mbh = lO^M© (see Fig- 
ure[2]). It can be found that our theoretical model also pre- 
dict that Rm is anti-correlate with i2-iokGv/iEdd as ob- 
servational features of LLAGNs. This is because the loud- 
ness Rm of our model scales Rm L-^/Lx oc m"^'^, while 
Lx/^Edd oc rh^. We find that the radio emission of most 
LLAGNs can be explained by our model if ^ 0.01 — 10% of 
the electron energy is ejected into the nonthermal electrons 
(see Figure [3]). However, pure RIAF without power-law 
electrons {rj = 0, see Figure [3]) always underpredicts the 
radio emission by several orders of magnitude compared 
the observations. 




-7 -6 -5 -4 



log(i2-10 kev/i'Edd) 

Fig. 3. — The relation between the mass-corrected radio loud- 
ness, i?M = i?x/Afgg, and the Eddington-scaled X-ray luminosity, 
i2-l0kov/iEdd. where Hx = /^5GHz/i2-l0kcV- The solid points 
represent the LLAGNs in our sample, and the solid lines represent 
the model prediction with the parameter r]= 10%, 1%, 0.1%, 0.01%, 
and 0% (from top to bottom) respectively. 



We further fit the radio to X-ray spectrum of all sources 
in our sample using the RIAF model with hybrid electrons. 
The results are shown in Figures HHHl and the adopted pa- 
rameters, rhout and rj, of each source are listed in Table 
(1). With exception of the optical/UV waveband in sev- 
eral sources, the model can successfully fit the SED of these 
sources from the radio to hard X-rays with only two free 
parameters. The X-ray emission dominantly come from 
the thermal electrons of the RIAF for all sources in our 
sample, while the observed radio emission is always 2-3 or- 
ders of magnitude higher than that predicted by the purely 
thermal RIAF model. The radio emission can be well de- 
scribed by the synchrotron radiation from the nonthermal 
electrons in the RIAF if a small fraction of the electron 
thermal energy resides in the power-law tail, where the rj 
range from 0.03% to 15%, with a typical value of ^^3%, for 
the sources in our present sample (see Table 1). 

5. SUMMARY AND DISCUSSION 

In this work, we explore the hybrid thermal-nonthermal 
synchrotron emission from the hot accretion flows of the 
RIAF. Such hot accretion flows are expected to exist in 
sources with low-mass accretion rates (e.g., LLAGNs). We 
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Fig. 4. — Models for the SED of LLAGNs. The solid points represent the observed emission, and the empty eircles denote the lower limit 
for the intrinsie AGN optical emission that derived from the variability of optical emission. The short-dashed and long-dashed lines represent 
the emission from the thermal and nonthermal electrons of the RIAF respectively, while the solid line show the radiation from both of these 
two components. 
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Fig. 5. — The same as Figure]?] but for the remaining 10 LLAGNs. 
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calculate the spectrum of the global RIAF model with hy- 
brid electrons, and find that synchrotron emission from 
the nonthermal electrons can be up to 2-3 orders of mag- 
nitude greater than that from the purely thermal electrons 
even with only several percent of electron thermal energy 
injected into the nonthermal electrons. This model has 
expl ained the multi-w avelength spectrum of Sgr A* very 
well (l Yuan et al.ll2003f ), and we further extend this model 
to other nearby LLAGNs, which have only a compact ra- 
dio core and have no evidence of the jet. Our model can 
roughly reproduce the observed anticorrelation between 
the mass-corrected radio loudness and Eddington ratio as 
found in the LLAGNs. We further perform the detailed 
modeling of the SED for 20 LLAGNs in our sample from 
the radio to X-ray waveband. We find that the SED of 
LLAGNs can be well described by our model, where the 
X-ray emission is produced predominantly by the inverse 
Compton scattering of the seed synchrotron photons pro- 
duced by both thermal and nonthermal electrons in the 
RIAF, while the radio emission mainly come from the non- 
thermal electrons. 

It is now established that RQ AGNs are not radio- 
silent, and do emit the radio emission at some level. 
iHo fc UlvestadI ()200lh found that 85% of the nearby 
Seyfert nuclei show the nuclear radio emission at 5 GHz, 
and their typical radio morphology is a compact core 
(eithe r unresolved or slightly resolved). lAnderson et al.l 
(|2004[ ) observed six nearby LLAGNs with high-resolution 
VLBA, and found that the radio emission is still unre- 
solved even at milliarcsecond scale, which roughly cor- 
responds to several to ten thousands gravitational radii. 
Therefore, contrary to the RL AGNs which have been ob- 
served with jet or core-jet structures directly, the phys- 
ical origin of the radio emission in RQ AGNs is much 
unclear. The most popular model is that there may ex- 
ist a scaled-down version of the jet in RQ AGNs. The 
coupled RIAF-jet model has been explore d to fit the mul- 



tiwavelength spectrum of LLAGN s (e.g., Wu et al] 2007t 



lYuan et al.l l2009al: lYu et all [2QTTI: iNemmen et al.l (20121) . 

We note that it is reasonable to apply the RIAF-jet model 
to the sources with observed jet structures (e.g., FR Is), 
while it is still controversial whether the same model can 
be used to the LLAGNs with only compact radio cores, 
where these sources have no evident jet structure even 
at scale of several thousand gravitational radii. Another 
possibility is that the radio emission in RQ AGNs also 
originate in the hot plasma (e.g., RIAF in LLAGNs and 
disk-corona in bright AGNs), where the turbulence, weak 
shocks and/or magnetic reconnection events may occasion- 
ally accelerate a fraction of the electrons to a harder power- 
law tail. The synchrotron emission from these nonthermal 
power-law electrons may account for the radio emission 
from the compact cores of these LLAGNs. The hybrid 
thermal-nonthermal synchrotron emission from the RIAF 
have well reproduced the two components of the radio 
spectrum and also the X-ray flares observed in Sgr A* (e.g., 
lYuan et al.l l2003h . We confirm the result that the pure 
RIAF alw ays underpredicts the radio emission for nearby 
LLAGNs (|Wu fc Caoll2"005h . However, the low-frequency 
radio emission of the compact cores in LLAGNs can be ac- 
counted if only a small fraction of the viscous dissipation 
energy in the accretion flow goes into accelerating electrons 



to a nonthermal power-law distribution (see Figures . 
It is well established that magnetic reconnection should be 
unavoidable and play an important role in co nverting the 
magnetic energy into the particles ( e.g., iHawie v & B albusI 
I2002t iGoodman fc Uzdenskvl I2008D Bo th the magnetic 
reconnection itself fe.g.. lDing et al.ll2010[ ) or the diffusive 
shock caused by the violent pl asma motion s in the mag- 
netic reconnection region (e.g.. ISpruitlll988D can naturally 
accelerate a small fraction of the electrons to be the power- 
law distribution. A small fraction of the power-law elec- 
trons presented in the RIAF will not affect the global 
structure of the RIAF (e.g., lOzel et al.|[2000[) . There- 
fore, our work provide the possibility that the radio emis- 
sion of nearby LLAGNs may originate in the nonther- 
mal electrons of the hot accretion flow. It may be the 
similar case for other bright RQ AGNs, where the radio 
emissio n is dominated by the nonthermal electrons in the 
corona. iLaor fc Behaij (j2008l) found the similarities of the 
radio/hard X-ray correlation in RQ AGNs and magneti- 
cally active stars, which support that the radio emission of 
RQ AGNs also dominantly come from the corona as that of 
active stars. We will further test this issue in a subsequent 
work for bright RQ AGNs considering the possible physical 
mechanism for the production of thermal and nonthermal 
el ectr o ns in the corona above the SSD. 

iHol (|2002f) found a strong anticorrelation between the 
radi o loudness and Eddington ratio for the LLAGNs (see 
alsolSikora et al.ll2007l : iPanessa et "all 120071: lYounes et all 
12012 ) . We investigate the dependence of the mass- 
corrected radio loudness, Rm = Rx/M^-^, on the Edding- 
ton ratio, i2-iokcv/iEdd, where the scaling of is de- 
rived from our model as that Lx oc Mbh and Lr oc Mg^. 
We find that the mass-corrected radio loudness is still anti- 
correlate with the Eddington ratio, where both quantities 
are not affected by the mass term. It is interesting to note 
that the dependence of the radio emission on the BH mass, 
Lr oc M^4, in onr model is similar to that of jet model 
(e.g.. iHeinz fc SunvaevI [2003h and that der ived from the 
observed fundamental plane of BH activity ([Merloni et al.l 
[2003; Falcke et al.,,2004i ). Therefore, it should be cautious 
to get the conclusion that the origin of the radio and X-ray 
emission is similar for the stellar-mass BH X-ray binaries 
(XRBs) and the supermassive BH AGNs based only on the 
fundamental plane equation of the BH activity, since that 
the dependence of the radio emission on the BH mass is 
similar for our model and the jet model. We find that the 
relation between the radio luminosity and X-ray luminos- 
ity is oc at a given BH mass in our model based on 
Lr oc rnP-^ and Lx oc Our model prediction on the 

radio— X-ray correlat ion is much shallower than that ob- 
served in XRBs re.g.. lGallo et al.ll2003l: [Corbel et al.ll2003D 
or that predicted by accretion-jet model (e.g.. lYuan fc Cuil 
;2005), where Lr oc Lx''- O^'" results provide a diagnostic 
that can distinguish the possible origin of the radio emis- 
sion in LLAGNs by investigating the radio— X-ray corre- 
lation for LLAGNs with the same or similar BH masses. 
However, our present sample is still limited, which prevent 
us to further test this issue. 

From a theoretical perspective, there are still some un- 
certainties in the current RIAF model with hybrid ther- 
mal and nonthermal electrons. The dependence of the 
spectrum on the model parameters have been explored in 
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forme r works (e.g., iQuataert fc Naravanlll999t lOzel et all 
l2000t) . The RIAF spectrum is not sensitive to some pa- 
rameters (e.g., a, /3), we adopt the typical values a = 0.3 
and /3 = 0.9, which are constrained from observations 
and/or simulations. Wind parameter s = 0.3 is adopted 
directly, which is constrained from observations of Sgr A*. 
We adopt the radius i?out = 5000i?g as the outer bound- 
ary. We find that the X-ray emission mainly come from 
the innermost region of the accretion flow, which is less 
affected by the outer boundary condition if it is larger 
than several tens gravitational radii. However, we find 
that the low-frequency radio emission come from the much 
larger regions from several tens to several thousand grav- 
itational radii. The LLAGNs in our sample have very low 
Eddington ratios (10^'' < ^2-10 kcv/-^Edd 5, 10"'*), and 
the whole accretion may be through RIAF. iLiu &: TaamI 
12009) proposed that the truncation radius should be 
around one thousand to several tens of thousands grav- 
itational radii (if the possible outer SSD is present) for 
m ~ 10^^ — 10^^ as adopted in fitting the SED of the 
LLAGNs in our sample. Therefore, the assumption of the 
outer boundary i?out = 5000i?g should be reasonable, and 
will not affect our main results. The region of the ra- 
dio emission from the nonthermal electrons in RIAF is 
more or less consistent with the size of the compact radio 
cores as observed b y the VLBI for some nearby LLAGNs 
(e.g.. IWrobel &: H0H2OO6;) . We neglect the contribution of 
the emission from the possible outer SSD (if present) in 
our SED fittings, since that we still don't know whether 
the outer SSD still exists in the LLAGNs with extreme 
low Eddington ratios or not, and the optical emission 
also easily suffer the contamination from the host galax- 
ies or contaminate by the nuclear starbursts. For exam- 
ple, IStorchi-Bergmann et al.l ()2005[ ) reported an evidence 
of recent starburst formation in the nucleus of NGC 1097, 
and its optical/UV emission may dominated by the stel- 
lar contribution. We find that our model prediction on 
the optical emission is slightly smaller than the observed 
emission but higher than the lower limits of the possible 
intrinsic AGN optical emission that derived from optical 
variabilities. Therefore, the observed optical emission of 



LLAGNs is roughly consistent with our model predictions 
even without considering the possible contribution from 
the outer SSD (e.g.. Figures |1][5]) . The parameter S, de- 
scribing the fraction of the turbulent dissipation that di- 
rectly heats the electrons in the flow, is still unclear, which 
may be in the range of several percent up to several tens 
perce nt based on recent simulations (e. g., Sharma et all 
'2007) and/or model fitting of the Sgr A* (lYuan et al.ll2003l 
i2006. ). We find that most of LLAGNs in our sample can 
be better fitted if (5 = 0.1 is adopted (see Figures |1][5]) . It 
should be noted that the adopted 5 = 0.1 i s slightly smaller 
than that used in Sgr A* {6 = 0.3. 1 Yuan et al. 2003), which 
may caused by different microphysics since that the accre- 
tion rate in Sgr A* much smaller than that used in our 
work. Our results will be roughly unc hanged even if w e 
use S = 0.3 as that used in Sgr A* (e.g.. lYuan et al.ll2006[ ). 
but the model fitting of the X-ray emission is less perfect 
as that with S — 0.1. The values of p and r] that describ- 
ing the nonthermal electrons is not strongly constrained, 
but which are degener ate for low-frequency radio emission 
(e.g.. lOzel et al.ll20Q0[ ). We fix the typical value of p=3.0, 
and then adjust the parameter rj to fit the radio emission of 
LLAGNs in our sample. We find that only a small fraction 
of the energy in nonthermal electrons is sufficient to pro- 
duce the radio emission at low frequencies. The value of 77 
can even be smaller if the parameter p is also smaller. Our 
fitting results show that the value of 77 in these LLAGNs 
of our sample is more or less co nsistent with that con- 
strained for Sgr A* (e.g., ~1% in lYuan et al"1l2003ll2006D 
and tha t derived from th e gamma-ray background (e.g., 
~4% in llnoue et al.l[2008l) . 
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